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Electronic, magnetic, and dynamical properties of the new superconducting /3-pyroclilore KOS2O6 
and related RbOs206 and CsOs2 06 compounds are calculated and compared with experiment and 
contrasted with structurally related spinel pyrochlores. The calculated susceptibility Stoner en- 
hancement (110%) and thermal mass enhancement A — 2.5-3 reflect moderate but perhaps impor- 
tant Coulomb correlations. The K"*" ion optic mode is found to be unstable, allowing large excursions 
of 0.5-0.6 A from its ideal site of the K ion along (111) directions. This dynamical mode is much less 
anharmonic in the isostructural Rb and Cs compounds (with larger cations), perhaps accounting 
for their progressively lower values of T^. Electron scattering from this very anharmonic mode may 
be the cause of the anomalous concave-downward resistivity that is seen only in KOsgOg. 

PACS numbers: 74.70.Dd,74.25. Jb,74.25.Kc 



I. INTRODUCTION 

Transition metal oxides with metal ions (T) lying on 
a pyrochlore sublattice display a wide variety of behav- 
ior. In the 1970s LiTi2 04 was the anomalous "high 
Tc" oxide superconductor amongst intermetallics, with 
its Tc ~ 13K. The vanadate LiV2 04 became the first 
(and still essentially the only) true heavy fermion metal 
based on d electrons rather than / electrons. Consider- 
able theoretical speculation on the microscopic basis for 
this heavy fermion behavior has left no consensus. Other 
members of the pyrochlore structure, such as CuIr2S4 
and TI2RU2O7, exhibit charge ordering and accompany- 
ing structural adjustment, possibly associated with the 
high symmetry of the undistorted structure. The "py- 
rochlore" sublattice occupied by the transition metal ions 
is comprised of corner-sharing tetrahedra that are known 
to be highly frustrating for nearest neighbor antiferro- 
magnetic (AFM) spin couphngs. 0, 

Recently Yonezawa et al.^ 01 reported the discov- 
ery of a superconductor KOS2O6 with a Tc = 9.6K, 
with /3-pyrochlore structure, a new variant of the py- 
rochlore class, and superconductivity was quickly ob- 
tained in isostructural and isoelectronic RbOs206 (6.3 
K)|| and CSOS2O6 (3.2 K)[| as well, suggesting new 
physics generic to this structural variant. The other ex- 
ample of superconductivity in Ad/hd pyrochlore oxides 
is Cd2Re207 with Tc = 1 K.Q There is no obvious in- 
dication of strong enhancement by correlation effects in 
the quasiparticle mass for the /3-pyrochlore compounds. 
Using the specific heat jump AC at Tc and the weak 
couphng relation AC/Tc = 1.437, Hiroi et al. 3 ob- 
tained a linear specific heat coefficient 7 = 19 mJ/K^ 
mol-Os, which is not especially large. The magnitude has 
been confirmed for RbOs206 by Briihwiler et al. who 
obtained 7= 17 mJ/K^ mol-Rb from the heat capacity. 
However, one highly unusual feature is that, while the re- 



sistivity p(T) of both RbOs206 (Ref. Q) and CsOsaOg 
(Ref. jg) have the normal upward curvature at low tem- 
perature, KOS2O6 with its higher Tc has a very peculiar 
concave downward shape'?] of p{T) immediately above 
Tc extending to 200 K. 

Conventional pyrochlore oxides have the chemical for- 
mula A2T2O7 or more descriptively A2T2O6O', where A 
is a larger cation, T is a smaller transition metal cation, 
and O' is an oxygen ion in a large tetrahedral site rather 
than forming the octahedron surrounding the transition 
metal ion. The distinguishing feature is the T pyrochlore 
sublattice, a network of corner-sharing tetrahedra that 
has been widely studied in the context of frustrated anti- 
ferromagnetism. Yonezawa et aZ.Q produced the (3 vari- 
ant with the general formula AT2O6 where A is a large 
monovalent alkaline metal cation. KOS2O6 has the same 
space group FdSm as conventional pyrochlore, but the 
large cations K+, Rb"*", Cs"*" are located (unexpectedly) 
at the O' site of the conventional pyrochlore structure, 
thus becoming the □2OS2O6K variant of conventional py- 
rochlore (where □ denotes an empty A site). KOS2O6 has 
the lattice constant 10.101 A and two formula units per 
fee primitive cell. 

Although the T ions lie on the same pyrochlore sublat- 
tice as in the closely related spinel system AT2O4, and are 
the centers of TOe octahedra, there are qualitative dif- 
ferences: in the spinel structure all the edge-sharing TOg 
octahedra are aligned along the cubic axes with T atoms 
being bridged by two O ions, while in the conventional 
pyrochlore and its /3 variant there are four different ori- 
entations of the vertex-sharing TOg octahedra and each 
pair of T ions is bridged by a single O ion. The four dif- 
ferent TOg octahedra are rhombohedrally distorted along 
the one of the four cubic body diagonals by an amount 
determined by the O internal parameter (u=0.3125 cor- 
responds to an equilateral octahedron). The pyrochlore 
lattice of the T sites can be viewed as a 3D generalization 
of the 2D Kagome lattice (see Fig. 



2 




FIG. 1: (color online) The crystal structure of the /3- 
pyrochlore superconductor KOS2O6. The largest atoms are 
K, the smallest dark spheres are O and are connected by 
sticks with the midsize atoms, Os. The pyrochlore sublattice 
formed by Os tetrahedra and truncated-tetrahedral cavities 
containing K ions is highlighted. 

In this paper we initiate the theoretical study of the 
KOS2O6 class of compounds, focusing on the electronic 
structure in comparison to the well studied conven- 
tional pyrochlore and spinel systems. The tight-binding 
description differs in important ways from the spinel 
(LiV2 04, say) counterpart, and the density of states 
N(Ei?) at the Fermi level Ep is not high. In fact it is 
below the average for the t2g band, due to Ej? lying in a 
valley, suggesting either electron or hole doping could en- 
hance Tc by raising N(Ei?). We further describe the dis- 
covery of an extreme dynamical instability corresponding 
to displacements of 0.5 A or more of the alkali ion along 
(111) direction. 

II. METHOD 

We have used both the full-potential nonorthogonal 
local-orbital minimum-basis (FPLO)scheme within the 
local density approximation (LDA)"^ and the full po- 
tential linearized augmented plane wave (LAPW) plus 
local orbitals method as implemented in Wien2k. [lol| 
The exchange and correlation potential of Perdew and 
Wang 11] was used. In FPLO, K 3s, 3p, 4s, 4p, 3d states, 
Os 4s, 4p, 4d, 4/, 5s, 5p, 6s, 6p and O 2s,2p, 3d were in- 
cluded as valence states. The inclusion of the relatively 
extended 3s, 3p, 4s, 4p, 3d semicore states as band states 
was done because of the considerable overlap of these 
states on nearest neighbors. The LAPW basis set is 



characterized by the atomic sphere radii of 2.05 bohr 
for K and Os and 1.4 bohr for O, plane-wave cut-off 
R-mtK„ia2;=5.5 and K-s, -p, Os-s, -p, -f and, 0-s local 
orbitals. The self-consistent potentials were calculated 
using 256 k points in the irreducible zone. The K, Os and 
O atoms are located at positions (|, |, |), (0, 0, 0),and 
(m, |, |) with site symmetries 43m, 3m, and mm2 respec- 
tively. We obtained the internal coordinate u — 0.317 by 
minimizing the total energy with FLAPW method. The 
corresponding deviation of the 0-Os-O bond angle from 
the right angle is ± 1.8°. Similar u values of 0.316 and 
0.315 were obtained for RbOs206 and CSOS2O6 respec- 
tively, corresponding to 0-Os-O bond angle deviations 
for right angle of 1.6° and 1.4°. Due to the relaxation of 
the O-Os-0 bond angle the relative differences of Os-O 
bond lengths are less the relative differences of the lattice 
constants for the three systems. 



III. RESULTS AND DISCUSSION 

A. Band and Density of States 

The guideline K+0s2^0g~ leads to the formal valence 
p = 5.5, or a d^'^ occupation (out of the six t2g states, 
including spin). As we see in Fig. |21 this puts E^? within 
the t2g bands, somewhat below half-filling of the threefold 
orbitally degenerate complex of bands. Surrounding Ep 
we have the t2g manifold (as in both pyrochlore A2T2O7 
and spinel AT2O4 systems) of width 3.2 eV (2.9 eV), 
and above it and separated by 1.6 eV (1.8 eV) is the 
eg manifold. The "O 2p bands" (with both Cg and t2g 
character mixed in) are fully occupied and separated from 
the t2g bands by 0.7 eV (1.5 eV) as shown in Fig. El 
(scalar- relativistic values are shown in the brackets). The 
non-cubic part of the 3m crystal field on the Os site is 
weak and the band structure shows no apparent further 
splitting of the t2g complex. We show a blowup of the 
t2g bands around Ep of KOS2O6 in Fig. [3 The twelve 
bands (3 t2g x 4 Os ions) are 5/12 filled. For comparison 
we have also calculated the band structures of the other 
/3-pyrochlore superconductors, CSOS2O6 and RbOs206. 
The band structures are extremely similar, with the same 
bandwidth, differing only in fine details. 

In Fig. 13 we display the density of states (DOS) of 
KOS2O6. The calculated Fermi level Ep lies within a 
valley, giving a value of N(Ei?)=8.3 states/eV/unit cell 
(scalar-relativistic value) that is slightly above the aver- 
age value of 7.5 states/eV in the t2g bands, and Ep lies 
0.3 eV above the peak in N(E) arising from the rather 
flat bands at the bottom of the complex. The DOS plot 
also indicates that the K states appear nowhere close to 
the Fermi level and do not exhibit mixing with the Os d 
states, which limits the role of the alkali ion to a donor 
of an electron and a source of electrostatic potential (and 
thus a possible scatterer). 
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FIG. 2: The Os t2g band complex of the superconductor 
KOS2O6 along fee symmetry lines. The formal d^'^ configu- 
ration results in 5/12 filling of these twelve bands. The rel- 
ativistic bands are marked with full (black) lines, the scalar- 
relativistic with dotted (red) line. 



B. Tight Binding Representation 

In order to understand the band structure of the t2g 
complex, we are thus left to consider the Os lattice 
bridged by O ions, suggesting a simple tight-binding de- 
scription. Crystal field splitting reduces the number of 
relevant orbitals to three t2g states per Os site. Sep- 
aration of the O bands suggests that the O p states 
can be "integrated out" and taken into consideration im- 
plicitly through the effective Os-Os hopping amplitude. 
Fujimoto l2] has shown that a single (s) state on a py- 
rochlore lattice with near neighbor coupling only leads 
to two dispersing bands and two absolutely flat bands, 
which reflects the close relation to the Kagome lattice 
(which has a single flat band). Singh et al. 13] consid- 
ered a much more realistic dda, ddn 12-band model for 
LiV204, and noted that for ddn =^dda, two of the twelve 
bands become flat. Unlike for the spinel structure, in the 
/3-pyrochlore structure the principle axes of the states 
are not aligned with the cubic axes so such a model does 
not apply directly to KOS2O6. 

Based on the geometry of the Os-O-Os bonds we have 
developed a tight-binding model with three t2g states on 
four Os sites in the unit cell. In an equilateral octahe- 
dron only dpi: hopping is possible for t2g orbitals (e.g. 
for O atom at (001) vertex of the local octahedron the 
dxz — > Px and dyz — » Py hopping is allowed by symmetry). 
The effective hopping amplitude between t2g orbitals on 
neighboring Os sites is then determined by this single pa- 
rameter and the relative orientation of the corresponding 
OsOg octahedra. This model (see Fig. 0} provides a rea- 
sonable band picture except that (1) the order of the T2g 
and Tig triplets (second and third levels from the bottom 
of the t2g complex) at the zone center are interchanged 
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FIG. 3: (color online) Scalar-relativistic density of states 
(solid line) and the site and orbital projected densities of 
states of KOS2O6: Os-t2g (dark shading), Os-eg (light shad- 
ing), O (dotted line), and K (dot-dashed line). In the inset we 
show in detail the effect of spin-orbit coupling to the density 
of states of the t2g manifold. 



and thus the band topology in this region (at and below 
Ei;-) is incorrect, (2) the four lower bands are too flat, and 
(3) unphysical degeneracies occur at the zone boundary 
points X and L (this last point is minor). Inclusion of 
the coupling to Cg states readily corrects the ordering of 
the Tig and T2g triplets and thus the band topology as 
shown in Fig. ^ (right panel) . Explicit incorporation of 
the O 2p states into a total 56 band model gives no ap- 
preciable change of the Os derived bands, indicating that 
the O 2p bands can be integrated out fairly accurately. 
The additional dispersion of the flat bands is likely to 
be due to direct Os-Os hopping. A rough estimate of 
this effect can be obtained from the Os lattice without 
O, yielding the d bandwidth less then 1 eV. 



C. Fermiology 

An appreciable effect of spin-orbit coupling is to be ex- 
pected in Os 5d bands. While the relativistic bands are 
qualitatively similar to the scalar-relativistic ones, there 
are some important differences. Removal of some band 
crossings leads to flatter bands and additional peaks in 
the density of states. As a result the density of states 
at the Fermi level is enhanced (see below) and the band 
structure in the vicinity of Fermi level becomes more sen- 
sitive to the oxygen position as we discuss below. 

There are two bands crossing the Fermi level in Fig. 
121 The flrst band gives rise to two closed sheets of the 
Fermi surface centered at the T point, which enclose the 
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FIG. 4: Tight-binding bandstructure: t2g-to-t2g hopping via 
oxygen only (left panel), including e^-to-eg and eg-to-t2g hop- 
ping (right panel). 




FIG. 5: The connected the Fermi surface with necks along 
the three- fold axis (spin-orbit coupling included). 



electron part of the reciprocal space in between of them 
(see Fig. EJ. The second band leads to a connected 
Fermi surface with necks along the three-fold axis (see 
Fig. |BJ. Visual inspection of the closed sheets of the 
Fermi surface suggests a possibility of partial nesting, in 
particular along the T — K direction. In order to pursue 
this possibility we have calculated the imaginary part of 
generalized susceptibility in the low energy limit 

Im x(q, Lo)^7rLuY, Sie(k) - (e(k + q) - Ep) (1) 

k 

= 7rcL'i/(q). (2) 

We emphasize that a large number of k-points is required 
to get reliable results, see the caption to Fig. [3 The cal- 
culated Imx(q) along T-K has some sharp Fermi surface 
related structure, but the peaks are only ±15% from a 




FIG. 6: The closed sheets of the Fermi surface centered on 
the r point (spin-orbit coupling included). 



smooth background and so do not suggest electronic in- 
stabilities. 

For the three /3-pyrochlores reported so far (K, Rb, 
Cs) IM IS the transition temperature Tj, and the lat- 
tice parameter vary inversely as shown in Table It is 
plausible to assume that the structural change is due to 
the size effect of alkaline metal ions, because the bands 
near the Fermi level consist of Os 5d orbitals with minor 
contribution from O 2p orbitals. The decrease of Tc with 
increasing volume under negative chemical pressure is in 
contrast to the case of conventional BCS superconductiv- 
ity in a single band model, where the Tc increases under 
negative pressure, because the density of states increases 
(see Table ID . 

Hiroi et al.^] have inferred the linear specific heat co- 
efficient for KOS2O6 of 7=19 mJ/K^ mole-Os from the 
heat capacity jump at Tc assuming the weak coupling 
relation AC/Tc = I.437. A similar value has recently 
been obtained by Briihwilcr et al. [i| for RbOs206 7=17 
mJ/K^ molc-Os directly from the heat capacity. The cal- 
culated value of N{Ef)=2S.2 states/Ry/Os for KOS2O6 
corresponds to a bare value 7(,=4.9 mJ/K^ mole-Os. The 
corresponding thermal mass enhancement A given by 



TABLE I: The experimental lattice constants and transition 
temperatures |^ 1^ together with the calculated densities 
of states at the Fermi level N(Ef) normalized per Os atom 
and the corresponding unrenormalized linear specific heat co- 
efficient 7i, per mole Os. 



Compound 


a(A) 




N(Ef)(Ry-^) 


76(mJ/K") 




KOsaOs 


10.101 


9.6 


28.2 


4.9 




RbOsaOs 


10.114 


6.3 


28.6 


5.0 




CsOsaOs 


10.149 


3.2 


30.3 


5.3 
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FIG. 7: The low frequency limit of the generalized sus- 
ceptibility i^{q) calculated along the F-K line obtained 
with 440x440x440 k-points (circles) in the cubic Brillouin 
zone. In the inset we compare with a curve obtained with 
160x160x160 k-mesh (diamonds). 

7/7fc 1 + A gives a large renormalization A — 2.9. The 
value A = 2.4 is obtained from the data for RbOs206. 
This renormalization includes the phonon contribution 
but most likely is due primarily to electronic processes. 

Important insight into magnetic fluctuation effects is 
provided by the enhancement of the bare Pauli suscepti- 
bility in a metal. We have evaluated, within the den- 
sity functional theory formalism, the Stoner enhance- 
ment of the susceptibility x = [i^in°{Ef)) = ' "^ti^re 
Xo = 'i^J^B^i^p) is the non-interacting susceptibility and 
S gives the electron-electron enhancement in terms of 
the Stoner constant /. We have calculated / using both 
the Janak-Vosko-Perdew theory and fixed spin mo- 
ment calculations, obtaining a value of S" = 2.15 ± 0.05, 
which does not indicate any ferromagnetic instability of 
the paramagnetic groundstate. The spin-orbit coupling 
was neglected in this calculation. 

D. Dynamical Instability 

In light of interesting but perhaps limited role of corre- 
lation effects (nothing like heavy fermion behavior), we 
have begun to pursue the character of electron-phonon 
coupling as a pairing mechanism. Since K"*" is a bare 
charge in a substantial hole in the /3-pyrochlore lattice, 
we have calculated the energy surface and deformation 
potential for a K-K "bond stretching" motion of the K 
ions, which lie on a diamond sublattice (of course, there 
is no K-K bond, the ions being ionized and also sepa- 
rated by cLk-k — -^o, = 4.33 A). The result wc find is 
a dynamical instability of the K ion. Although the force 
vanishes by symmetry for the ideal structure, for increas- 
ing separation of K ions lying along the (111) direction 
the energy decreases. The energy is minimized only after 
a displacement of the K ion by 0.65 Al This motion is 
directed along (111) channels in the OS2O6 system, and 



this crystal structure may not be stable for the smaller 
alkali cations Na and Li simply because they do not stay 
put near the ideal site. In Fig. |Hlwe show the energy as 
a function of the alkali ion displacement in this mode for 
K, Rb and Ce as well as fictitious Na compounds (the lat- 
tice constant of KOS2O6 was used for the Na compound). 
Note that the curves are upper bounds, since allowing the 
Os and O atoms to relax at any displacement would only 
lower the energy. 

The high symmetry position is found to be unstable for 
Na, while a very flat energy surface is found for K over 
a large range of displacements. The Rb system exhibits 
significant anharmonicity, which is reduced when going to 
Cs. Large differences in the energy surface between sys- 
tems which have very similar lattice constants and band 
structures can be understood as follows. If only the effect 
of the electrostatic potential on the nucleus at the alkali 
metal site was considered (i.e. if the electron charge was 
frozen) the site would be dynamically unstable although 
the force vanishes. This is reflected in the fact that the 
first non-spherical term in the site expansion of the elec- 
trostatic potential is a cubic polynomial, indicating an 
inflection point. The site is thus stabilized due to elec- 
tronic relaxation which is accomplished by mixing of the 
outer alkali ions orbitals with the orbitals on its neigh- 
bors. This explains the pronounced difference between 
Na, for which the 2p orbitals are substantially more lo- 
calized than the 5p orbitals of Cs. Note that shape of 
the instability or anharmonicity corresponds with the ef- 
fect of the electrostatic repulsion of the four neighboring 
alkali ions (which form a tetrahedron), which tends to 
move the atom in the center away from the vertices. 

In addition to the alkali ion displacement we have in- 
vestigated the symmetric O mode, which corresponds to 
varying the internal parameter around its equilibrium 
value of 0.317. This Raman-active mode correspond to a 
simultaneous rhombohedral distortion of the OsOg octa- 
hedra along a cubic body diagonal (different diagonal for 
each of the four Os atoms in the unit cell). Locally the O 
atom is moving perpendicular to the line connecting its 
nearest-neighbor Os pair. The calculated frequency is 65 
meV=525 cm~^. Similar frequencies of 65 meV and 64 
meV were obtained for RbOs2 06 and CSOS2O6 respec- 
tively. The deformation potential for the band crossing 
the Fermi level calculated at L point (very near the Fermi 
level) amounts to Aefc/Ai?=1.8 eVA"-'^, where AR is the 
displacement of each O ion. The presence of two bands 
in the vicinity of the Fermi level (at the L point and close 
to the center of F-L line) makes the Fermi surface rather 
sensitive to this oxygen mode. These two bands move 
in mutually opposite directions as the oxygen is moved 
away from equilibrium. As a consequence displacement 
of the O atom by less then 0.05 A results in appearance 
of an additional electron pocket centered at the L point 
followed by sticking together and opening of holes along 
F-L direction in the sheets centered at F point. 
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FIG. 8: (color online) Energy surface for a displacement A of 
alkali ions along (111) direction. The curves were obtained 
by spline interpolation of energies and forces evaluated at 
approximately ten points. Note the very large (~ 0.5 eV) 
instability of the Na ion. 

IV. DISCUSSION AND SUMMARY 

In this paper we have analyzed the electronic structure 
of the /?-pyrochlore KOS2O6, which is nearly identical to 
those of the Rb and Cs compounds. The Os t2g states are 
well separated from the more tightly bound O 2p states, 
and also separated from the unoccupied Cg states, leav- 
ing the focus on the complex of twelve t2g-derived bands. 
Spin-orbit coupling has a large effect on N(Ei;-), increas- 
ing it by 60% over its value if S-0 coupling is neglected 
(other regions will have decreased values of N(E)). This 
system has often been compared to the other pyrochlore 
structure superconductor Cd2Re2 07, which has Tc = 1 
K. The tig bands of the two compounds have very much 
the same shape, with those of KOS2O6 being 20% wider 
[TEll . The different band filling of these compounds pre- 
cludes serious comparison of their superconducting be- 
haviors, and the more interesting comparison is the fac- 
tor of three variation in Tc within the /3-pyrochlore class 
when they all have very similar electronic structures. 

A single parameter tight binding model provides a 
starting point for an effective model Hamiltonian, but 
leaves the lower four bands flatter than observed. This 
aspect, as well as ordering of bands at k=0, can be im- 
proved by including coupling between the t2g and Cg 
states. The calculated Fermi surfaces consist of two 
closed surfaces at the zone center, with flat portions, and 



a large surface that is open along the (111) directions. 
The closed surfaces give rise to nesting features at spe- 
cific values of q along the (110) direction, but while sharp 
the features are not very strong. 

Investigation of the strength of electronic correlations 
has been initiated here. Comparison of the measured 
linear specific heat coefficient 7 for RbOs206 with the 
calculated value of N(Ei?) leads to a dynamic quasipar- 
ticle mass enhancement A = 2.4, indicating important 
correlation effects but far from heavy fermion type of be- 
havior. Calculation of the Stoner enhancement of the 
q=0 susceptibility gives a factor of two, again indicating 
important correlation corrections but nothing close to a 
ferromagnetic instability. 

The only substantial difference between the K, Rb, 
and Cs compounds that we have found is the degree of 
(in)stability of the alkali cation in its tetrahedral inter- 
stitial site. The potential for K ion motion along the 
(111) directions is flat over a distance of almost 1 A, in- 
dicating an extremely floppy ion whose motion may pro- 
vide the scattering of carriers that is reflected in peculiar 
concave-downward resistivity at low temperature. The 
Rb and Cs ions are progressively more stable, and their 
resistivity behavior is more conventional. Calculation for 
the Na compound (not yet reported) indicates a seriously 
unstable potential surface, with the minimum being more 
than 1 A away from the normal high symmetry site of 
the cation. This instability may make the Na compounds 
unstable (and unsynthesizable) in this structure. 

The close similarity of the electronic structures of the 
compounds in the KOS2O6 series, as well as almost iden- 
tical equilibrium position and dynamics of the O ion, 
allows two explanations of the largely different Tc's: (i) 
fine details of the electronic structure at the Fermi level 
are very important, (ii) the very different degree of an- 
harmonicity of the alkali ion is responsible for the differ- 
ences in Tc. We propose that the pressure dependence 
of Tc can resolve these two scenarios. In scenario (i) a 
strong pressure dependence of Tc with a negative slope 
is expected following the trend accross the K-Rb-Cs se- 
ries. On the other hand the alkali ion dynamics, deter- 
mined primarily by its ionic radius, is not sensitive to 
small changes of the volume and so no dramatic pressure 
dependence of Tc is be expected if scenario (ii) applies. 
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01ER45876, while W. E. P. was supported by NSF Grant 
DMR-0421810. 
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